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, The ordering property of the classical pyrochlore Heisenberg antiferromagnet with the ferromag- 

■ netic next-nearest-neighbor interaction is investigated by means of a Monte Carlo simulation. The 

' model is found to exhibit a first-order transition at a finite temperature into a peculiar ordered state. 

While the spin structure factor, i.e., the thermal average of the squared Fourier amplitude of the 

. spin, exhibits a finite long-range order characterized by the commensurate spin order of the period 

' four, the thermal average of the spin itself almost vanishes. It means that, although the amplitude 

I of the spin Fourier component is long-range ordered, the associated phase degree of freedom remains 

, to be fluctuating. 

(N , 
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Recently, there arises a considerable interest in the ordering of pyrochlore magnets. The pyrochlore lattice consists 
of corner-sharing tetrahedra, and is generally regarded as a typical realization of geometrically frustrated lattices. In 
particular, the classical Heisenberg antiferromagnet on the pyrochlore lattice is known to exhibit no magnetic long- 
range order of any kind even at zero temperature, if the interaction is restricted only to nearest neighbors QIMH^Q- 
^ , This is due to the extremely high degeneracy of the ground state induced by severe geometrical frustration. Such a 
high degeneracy is realized on the delicate balance among frustrated interactions, and might be lifted by the weak 
perturbation which inevitably exists in real magnets, e.g., the further neighbor interactions, the dipolar interaction, 
the magnetic anisotropy, the quenched randomness, or the lattice distortion etc. The lifting of the degeneracy 
] would eventually lead to a phase transition at a low but nonzero temperature, possibly of an exotic type peculiar 
' to geometrically frustrated magnets. Thus, it is important to clarify the role of various weak perturbations on the 
O ' ordering of the nearest-neighbor Heisenberg model. 

Some time ago, Reimers performed a pioneering study of the ordering of the classical Heisenberg pyrochlore an- 
tiferromagnet with the antiferromagnetic first-neighbor interaction J and the further- neighbor interactions (up to 
the fourth- neighbors), J2 ~ t/4, within a mean-field approximation ^5]. He constructed a possible phase digram as a 
^ ' function of the interaction parameters J2/\J\ ~ J4/I J|, revealing a variety of ordered phases. 

' In the present paper, we investigate the ordering property of the classical pyrochlore Heisenberg antiferromagnet 
with the ferromagnetic next-nearest-neighbor interaction J2 > by means of a Monte Carlo simulation |^. In fact, 
Q,^ we also studied the same model with the antiferromagnetic next-nearest-neighbor interaction J2 < 0, to find that 
I the system exhibited a first-order transition into the q ~ ordered phase with the coUinear up-up-down-down spin 
\^ ' structure 0,0- There, the ordered state is conventional in the sense that the long-range magnetic order sets in. In 
the case of the ferromagnetic next-nearest-neighbor interaction J2 > 0, by contrast, while the system also exhibited a 
first-order transition at a finite temperature into the ordered state, this ordered state turns out to be quite peculiar, 
' as we shall see below. 



II. MODEL AND PHYSICAL QUANTITIES 



o 

O ' The model we consider is the classical Heisenberg model on a pyrochlore lattice with the antiferromagnetic first- 
neighbor interaction J < and the ferromagnetic second-neighbor interaction J2 > 0. The Hamiltonian is given 
by 

^: H = -j J2 s,-s,~j2 Sk-si, (1) 

<i,j> <k,l> 

where Si is a three-component classical Heisenberg spin variable at the i-th site. Si — {Sf, Sf , Sf) with \Si\ = 1, and 
the sum is taken over all first-neighbor and second-neighbor pairs on the lattice. 

We study the ordering properties of the model by Monte Carlo simulations. The strength of the second-neighbor 
interaction is set to J2 = 0.1| J|. (In fact, we have also studied other values of J2, to find that the result is insensitive 
to the J2-value, at least qualitatively.) At lower temperatures, a very slow dynamics, presumably originated from the 
intrinsic frustration effect, arises. To circumvent the resulting thermalization problem, we adopt in our simulation a 
temperature-exchange method combined with the standard heat-bath updating. 
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In a cubic unit cell of the pyrochlore lattice, there are 16 sites (spins). We measure the linear size of our system L 
in units of the cubic unit cell, i.e., our system contains TV = IQL^ spins in total. Periodic boundary conditions are 
imposed in all directions. 

On a pyrochlore lattice, complex spin order could arise due to the severe frustration effect. Accordingly, in order 
to detect the possible complex spin order, we calculate the spin freezing parameter, a quantity familiar in the study 
of spin glasses, defined by 

e^-Y.^q%), {p^.v = x,v,z\ = (2) 

where (• • •) is the thermal average, and the upper suffixes (a) and (b) denote the two copies (replicas) of the system. 
In order to compute gs , we simulate two independent replicas of the system, (a) and (b), by using different initial 

(2) 

conditions and different random-number sequences. In equilibrium, ql is identically given by 

'if'=w^i:^^^-s3?- (3) 

("2) 

As is evident from Eq.(3), the spin freezing parameter ql becomes unity for any ordered state where the mutual 
direction of spins is completely frozen, whatever type of ordered state it may be. If the mutual direction of spin 

(2) . (2) 

remains fluctuating, on the other hand, qs tends to vanish. Note that q^ is invariant under continuous symmetry 
operations of the Heisenberg Hamiltonian, i.e., under global spin rotations made independently of the two replicas (a) 

f 2) 

and (b). Such an invariance of ql under global spin rotations is practically important, since, in finite-size simulations 
of the isotropic Heisenberg model, global spin rotations easily occur in the course of simulation. Concerning the 

(2) 

discrete symmetry of the Hamiltonian, ql is invariant under global spin refiections, but is not necessarily so under 
discrete symmetry operations associated with the underlying lattice symmetry, e.g., lattice translations, which can 
not be adsorbed into spin rotations nor reflections. 

In order to monitor the type of spatial spin order, we also calculate the spin structure factor F{q), i.e., the thermal 
average of the squared Fourier amplitude defined by 

Pi^-^{\Y.^3e''-''\") (4) 
j 

Other standard quantities, e.g., the energy, the specific heat and the susceptibility, etc are also calculated. 



III. MONTE CARLO RESULTS 



In Fig.l, we show the temperature and size dependence of the energy per spin. As can be seen from the figure, an 
almost discontinuous jump characteristic of a first-order transition is observed at T = Tpf ~ 0.094J, which becomes 
increasingly eminent for larger lattices. The calculated energy histogram (not shown here) exhibits a pronounced 
double peak structure at around T — Tjv, demonstrating that the transition is certainly first order. 

In Fig. 2, we show the temperature and size dependence of the freezing parameter qs . As can be seen from the 
figure, the growth of gs^"* is rather week; gs^^ takes small values of order 10~^ even below the first-order transition 

f2) 

point. Indeed, even when we extrapolate the observed temperature dependence of ql down to T = 0, it yields only 

(2) 

10% of its full value or less. Whether qs in the ordered state entirely vanishes or keeps a small nonzero value in 
the thermodynamic limit L — > oo, is not entire clear. While the data of larger lattices L — 8 and 12 suggest that a 
small nonzero value persists even in the thermodynamic limit, further careful study is necessary to clarify this point. 
Anyway, the spin remains largely fiuctuating even in the ordered state. 

We also calculate the spin structure factor F{q) at wavevectors q = ^{h,h,l) at a temperature T = 0.084|J| 
below T/v (a being the spacing of the cubic unit cell), which is shown in Fig. 3 in the {l,h) plane. Main peaks are 
observed at (±5/4, ±5/4, 0) and at (±3/4, ±3/4, 2), with sub-main peaks at (±3/4, ±3/4, 0) and at (±5/4, ±5/4, 2). 
Note that, in the whole q-space, there are twelve independent points at each (5/4,5/4,0), (3/4,3/4,2), (3/4,3/4,0) 
and (5/4, 5/4, 2), respectively. These peaks turn out to be sharp, with very little amplitude left at other positions of 
q. The structure factor along the {h, h, 0) direction is shown in Fig. 4 for various lattice sizes L. These data clearly 
indicate that, when probed via the structure factor, the system exhibits an almost complete commensurate magnetic 
long-range order of the period four. We note that such a period-four commensurate spin order is at odd with the one 
predicted in Ref.(5| based on the mean-field analysis. 
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FIG. 1: The temperature and size dependence of the energy per spin. 
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FIG. 2: The temperature and size dependence of the spin freezing parameter q\ 
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Our observation for the spin structure factor is in apparent contradiction to the observed near-vanishing of the spin 

(2) 

freezing parameter qs ■ It means that, although the spin itself remains largely fluctuating even in the ordered state, 
the system is almost fully ordered when probed via the spin structure factor. Since the structure factor measures 
only the amplitude of the spin Fourier component irrespective of the state of its phase degree of freedom, the ordered 
state is a peculiar one where the Fourier amplitude of the spin exhibits a long-range order while the phase of the spin 
Fourier component remains fluctuating. 

We also examined the snapshot of spin patterns realized at lower temperature, to find that the period-four periodicity 

(2) 

of the ordered state is hardly visible in the spin pattern itself consistently with our observation of , which, however, 
might partly be due to the finite-temperature effect. 

(2) 

One may wonder if the observed near-vanishing of might still be an artifact due to finite-size effects. For 

(2) 

instance, although ql is invariant under both global spin rotations and reflections as mentioned above, possible 

(2) 

global sliding motion of the ordered state associated with lattice translations might tend qs to vanish. In order to 

(2) 

examine such a possibility, we also calculate qs by applying an appropriate symmetry-breaking field of intensity Hsb 
to the system, —HsbJ^i^i ' M^j where the direction of the symmetry-breaking field hf^ (with \hf^\ — 1) is chosen 
to be the spin direction in the lowest-temperature configuration generated in separate simulation with Hst = 0. If 
the observed near-vanishing of qs is a finite-size effect associated with the global sliding of the ordered state, the 
application of the symmetry-breaking field would quickly pin down such a global sliding, leading to the appearance 

(2) 

of larger qs values of order unity when the temperature is lowered across the transition temperature of the Hsb = 

(2) 

model. Our calculation with Hgb = 0.1| J2I, however, indicates that the induced qs still remains small, without any 

(2) 

rapid growth being observed around the transition temperature of the Hsb = model: The qs value extrapolated 
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FIG. 3: The spin structure factor F{q) at wavevectors q= ^{h, h, I) shown in the (/, h) plane. The temperature is T = 0.084| Jj 
below Tm- 
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FIG. 4: The /i-dependence of the spin structure factor F{q) along the {h, h, 0) direction for various lattice sizes L. The 
temperature is T = 0.084| J|. 



to T = is about 0.3 for L — A. These observations might suggest that the observed near vanishing of gs is a bulk 
phenomenon, not merely a finite-size effect due to the global sliding of the ordered state. 

IV. SUMMARY 

The ordering property of the classical pyrochlorc Heisenberg antiferromagnet with the ferromagnetic next-nearest- 
neighbor interaction is investigated by means of a Monte Carlo simulation. The model turns out to exhibit a first-order 
transition into a peculiar ordered phase where the spin structure factor, i.e., the thermal average of the squared Fourier 
amplitude, exhibits a finite long-range order characterized by the commensurate spin order of the period four, while 
the thermal average of the spin itself almost vanishes. 

We finally mention that an exotic ordered state as revealed here might have important implications to experiments. 
In such an exotic ordered state, while the conventional Bragg peak would be observed in neutron-scattering measure- 
ments, no magnetic long-range order would be observed in NMR or /iSR measurements. Experimental quest for such 
an exotic ordered state in real pyrochlorc materials would be of utmost interest. 

The author is thankful to Dr. M.J. P. Gingras for useful discussion. 
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